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Abstract 

We present the expansion of the two-pion exchange three-nucleon potential (TPE-3NP) to chiral 
order q 4 , which corresponds to a subset of all possibilities at this order and is based on the 7rN 
amplitude at 0(q 3 ). Results encompass both numerical corrections to strength coefficients of 
previous 0(q s ) terms and new structures in the profile functions. The former are typically smaller 
than 10% whereas the latter arise from either loop functions or non-local gradients acting on the 
wave function. The influence of the new TPE-3NP over static and scattering three-body observables 
has been assessed and found to be small, as expected from perturbative corrections. 

PACS numbers: 13.75. Cs, 21.30.Fe, 13.75. Gx, 12.39.Fe 
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I. INTRODUCTION 



The research programme for nuclear forces, outlined more than fifty years ago by Taketani, 
Nakamura, and Sasaki [1], treats pions and nucleons as basic degrees of freedom. This insight 
proved to be very fruitful. On the one hand, it implies the interconnection of all nuclear 
processes, both among themselves and with a class of free reactions. On the other, it 
determines a close relationship between the number of pions involved in a given interaction 
and its range. As a consequence, the outer components of nuclear forces are dominated by 
just a few basic subamplitudes, describing either single (N — > ttN) or multipion (nn — > tttt, 
ttN — > irN, ttN — > irirN, ...) interactions. 

Nevertheless, it took a long time for a theoretical tool to be available which allows the 
precise treatment of these amplitudes. Nowadays, owing to the development of chiral per- 
turbation theory (ChPT) in association with effective lagrangians [2, 3], the roles of pions 
and nucleons in nuclear forces can be described consistently. The rationale for this approach 
is that the quarks u and d, which have small masses, dominate low-energy interactions. 
One then works with a two-flavor version of QCD and treats their masses as perturbations 
in a chiral symmetric lagrangian. The systematic inclusion of quark mass contributions is 
performed by means of chiral perturbation theory, which incorporates low-energy features 
of QCD into the nuclear force problem. In performing perturbative expansions, one uses a 
typical scale q, set by either pion four-momenta or nucleon three-momenta, such that g< 1 
GeV. 

Nuclear forces are dominated by two-body (NN) interactions and leading contributions 
are due to the one-pion exchange potential (OPEP), which begins [4] at O(q ). The two- 
pion exchange potential (TPEP) begins at 0(q 2 ) and, at present, there are two independent 
expansions up to 0(q A ) in the literature, based on either heavy baryon [5] or covariant [6, 7] 
ChPT. The TPEP is closely related with the off-shell nN amplitude and, at this order, 
two-loop diagrams involving intermediate nir scattering already begin to contribute. 

In proper three-nucleon (3iV) interactions, the leading term is due to the process known 
as TPE-3NP, in which the pion emitted by a nucleon is scattered before being absorbed 
by another one. It has been available since long [8-10], involves only tree-level interactions 
and has the longest possible range. This contribution begins at 0(q 3 ) and consistency with 
available NN forces demands the extension of the chiral series for the 3NP up to 0(q 4 ). 
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However, the implementation of this programme is not straightforward, since it requires 
the evaluation of a rather large number of diagrams. With the purpose of exploring the 
magnitude of 0(q 4 ) effects, in this work we concentrate on the particular subset of processes 
which still belong to the TPE-3NP class. Our presentation is divided as follows. In section 
II we display the general relationship between the TPE-3NP and the irN amplitude, in order 
to discuss how it affects chiral power counting in the former. The irN amplitude relevant for 
the 0(q 4 ) potential is derived in section III and used to construct the three-body interaction 
in section IV. We concentrate on numerical changes induced into both potential parameters 
and observables in sections V and VI, whereas conclusions are presented in section VII. 
There are also four appendices, dealing with kinematics, 11N subthreshold coefficients, loop 
integrals and non-local terms. 

II. GENERAL FORMULATION 

Potentials to be used into non-relativistic equations can be derived from field the- 
ory by means of the T-matrix. In the case of three-nucleon potentials, one starts 
from the non-relativistic transition matrix describing the process N(pi) N{j>2) N(p$) — > 
N(pi) N(p' 2 ) N(p' 3 ), which includes both kernels and their iterations. The former corre- 
spond to proper interactions, represented by diagrams which cannot be split into two pieces 
by cutting positive-energy nucleon lines only, whereas the latter are automatically gener- 
ated by the dynamical equation. Therefore, just the kernels, denoted collectively by t 3 , are 
included into the potential. 

The transformation of a T-matrix into a potential depends on both the dynamical equa- 
tion adopted and conventions associated with off-shell effects. The latter were discussed 
in a comprehensive paper by Friar [11]. Here we use the kinematical variables defined in 
Appendix A and relate £3 to the momentum space potential operator W by writing [12] 



(P1.P2.P3 \W\p 1 ,p 2 ,p 3 ) 



(27i) 3 5 3 (P'-P) i 3 (Pi,P2,P' 3 ,Pi,P2,P 3 ) ■ 



(1) 



In configuration space, internal dynamics is described by the function 




xe 



i[Q r -(r'-r)+ Q p -(p'-p)+q r -(r>+r)/2+q p (p>+p)/2\ 




(2) 
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which is to be used in a non-local version of the Schrodinger equation: 



--v r 2 , 

m 



m 



\/3/2l drdpW(r',p';r,p)i/;(r,p) . (3) 



Non-local effects are associated with the variables Q r and Q p . When these effects are not too 
strong, they can be represented by gradients acting on the wave function and the potential 
W is rewritten as 
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W(r',p';r,p) =6 3 {r'-r)5 3 (p'-p) 2/VH V(r,p 



(4) 



The two-pion exchange three-nucleon potential is represented in Fig. la. It is closely related 
with the ttN scattering amplitude, which is 0(q) for free pions and becomes 0(q 2 ) within 
the three-nucleon system. As a consequence, the TPE-3NP begins at 0(q 3 ) and, at this 
order, it also receives contributions from interactions (c) and (d), which have shorter range. 
The extension of the chiral series to 0(q 4 ) requires both the inclusion of single loop effects 
into processes that already contribute at 0(q 3 ) and the evaluation of many new amplitudes, 
especially those associated with diagram (b). 



j. x_ 




(a) 



(b) 



(c) 




FIG. 1: (Color online) Classes of three-nucleon forces, where full and dashed lines represent nucleons 
and pions respectively; diagram (a) corresponds to the TPE-3NP. 

In this paper we concentrate on the particular set of processes which belong to the TPE- 
3NP class, represented by the T-matrix T nn and evaluated using the kinematical conditions 
given in Fig. 2. The coupling of a pion to nucleon % — (1, 2) is derived from the usual lowest 
order pseudo-vector lagrangian £W an d the Dirac equation yields the equivalent forms for 
the vertex 

(W 2 /*) V u (P'-P) 75 u] M = (mg A /f n ) [t u] {t) , (5) 
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where gA, fw and m represent, respectively, the axial nucleon decay, the pion decay and the 
nucleon mass. 




FIG. 2: (Color online) Two-pion exchange three-nucleon potential. 



The amplitude for the intermediate process 7i a (k)N(p) 
structure 

T ba = 5 a b T + + ie bac r c T~ 

and Fig. 2 yields 



n b (k')N(p') has the isospin 



(6) 



T = - 

■*■ 7T7T 



mg A 



[m7 5 m] (1) [n7 5 Ml (i 



k 2 — ji 2 k' 2 — fi 2 

[ T (l). T (2) T+ _ ZT (l) XT (2). T (3) T -](3) ? (7) 

fx being the pion mass. Results in Appendix A show that [m7 5 «]^ — > 0(q), whereas pion 
propagators are 0(q~ 2 ). As a consequence, in the 0(q A ) expansion of the potential one 
needs T nn to 0(q) and T to 0(q 3 ). For on-shell nucleons, the sub amplitudes T^ 1 can be 
written as 



HP 



D ± - —a^p'-pYK^B 



HP) 



(8) 



with K = (k'+k)/2. The dynamical content of the ttN interaction is carried by the functions 
D and B ± and their main properties were reviewed by Hohler [13]. The chiral structure 
of these sub amplitudes was discussed by Becher and Leutwyler [14, 15] a few years ago, 
in the framework of covariant perturbation theory, and here we employ their results. As 
far as power counting is concerned, in Appendix A one finds [u(p') u(p)]^ — > O(q ) and 
iiT^ip') cr^vip' —pYK u u(p)}^ — > 0(q 2 ), indicating that one needs the expansions of D ± 
and B ± up to 0(q 3 ) and 0(q) respectively. 



At low and intermediate energies, the 7rN amplitude is given by a nucleon pole superim- 
posed to a smooth background. One then distinguishes the pseudovector (PV) Born term 
from a remainder (R) and writes 

t ± = t± + t£. (9) 

The former contribution depends on just two observables, namely the nucleon mass m and 
the 7rN coupling constant g, as prescribed by the Ward-Takahashi identity [16]. The calcu- 
lation of these quantities in chiral perturbation theory may involve loops and other coupling 
constants but, at the end, results must be organized so as to reproduce the physical values 
of both m and g in T^ v [17]. For this reason, one uses the constant g, instead of (^a/Ztt), 
since the former is indeed the observable determined by the residue of the nucleon pole 
[13, 15, 18]. The pv Born sub amplitudes are given by 

g 2 ( k'-k k'-k 
2m V s — m 2 u — m 2 



s i = -» 2 — — — — . (ii) 

' \s — m z u — wr J 

D -.|L(J^_ J^_*), (12) 

' 2m \ s — m A u — m z m J 

B vv = -a 2 + + T-2) , (is) 

' \s — m z u — m z 2m z J 

where s and u are the usual ttN Mandelstam variables. In the case of free pions, their chiral 
orders are respectively [D+ v , B+ v , D~ v , B~ v ] — > 0[q 2 , g" 1 ,q, q°], but important changes do 
occur when the pions become off-shell. 

The amplitudes Tjjr receive contributions from both tree interactions and loops. The 
former can be read directly from the basic lagrangians and correspond to polynomials in 
t = (k' — k) 2 and v = (p'+p)-(k'+k)/Am, with coefficients given by renormalized LECs [15]. 
The latter are more complex and depend on Feynman integrals. In the description of ttN 
amplitudes below threshold, one approximates both types of contributions by polynomials 
and writes [13, 19] 

X R = J2xmnV 2m t\ (14) 

where Xr stands for D^, B^/v, D^jv or B^. The subthreshold coefficients x mn have the 
status of observables, since they can be obtained by means of dispersion relations applied 
to scattering data. As such, they constitute an important source of information about the 
values of the LECs to be used in effective lagrangians. 
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The isospin odd subthreshold coefficients include leading order terms, which implement 
the predictions made by Weinberg [20] and Tomozawa [21] for ttN scattering lengths, given 
by 

_ v 1 

Dy/T = ' WT = ■ (15) 

For free pions, one has [D^ T , B^ T ] — > 0[q, q°], but these orders of magnitude also change 
when pions become virtual. 

Quite generally, the ranges of nuclear interactions are determined by t-channel exchanges. 
At 0(q 3 ), the TPE-3NP involves only single-pion exchanges among different nucleons and 
has the longest possible range. Another t-channel structure becomes apparent at 0(q A ), 
associated with the pion cloud of the nucleon, which gives rise to both scalar and vector 
form factors [18]. These effects extend well beyond 1 fm [22, 23] and a limitation of the 
power series given by Eq. (14) is that they cannot accommodate these ranges, since Fourier 
transforms of polynomials yield only (^-functions and its derivatives. In the description of the 
irN amplitude produced by Becher and Leutwyler [15], one learns that the only sources of 
medium range (mr) effects are their diagrams k and /, which contain two pions propagating 
in the t-channel. In our derivation of the TPE-3NP, the loop content of these diagrams is not 
approximated by power series and, for free pions, the non-pole subamplitudes are written as 

d + r = D+ r (t) + [d+ + d+y + d+j] (2) + [d+y + d\ x vH + d+t 2 ] (3) , (i6) 

B+ = B+ r (t) + [b+u] {1) , (17) 
Dr = D mr (t) + [v/{2fD] (1) + [d m u + d w u 3 + d m vt] {3) , (18) 
B R = B~ r (t) + [1/(2/2) + g-j ^ + ^ b -y + g- ^ ^ ^ (19) 

where the labels (n) outside the brackets indicate the presence of 0(q n ) leading terms 
and mr denotes terms associated with the nucleon pion cloud. The bar symbol over some 
coefficients indicates that they do not include both Weinberg- Tomozawa and medium range 
contributions, which are accounted for explicitly. The functions D R and B R depend on the 
parameters /„., g^, /•*, m and on the LECs q and di, which appear into higher order terms 
of the effective lagrangian. The subthreshold coefficients are the door through which LECs 
enter our calculation and their explicit forms are given in Appendix B. 

The dynamical content of the 0(q 3 ) ttN amplitude is shown in Fig. 3. The first two 
diagrams correspond to PV Born amplitudes, whereas the third one represents the Weinberg- 
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FIG. 3: (Color online) Representation of the irN amplitude used in the construction of the TPE- 
3NP. 



Tomozawa contact interaction, all of them with physical masses and coupling constants. The 
fourth graph summarizes the terms within square brackets in Eqs. (16-19) and depends on 
the LECs. Finally, the last two diagrams describe medium range effects owing to the nucleon 
pion cloud, associated with scalar and vector form factors. This decomposition of the nN 
amplitude has also been used in our derivation of the two-pion exchange components of the 
NN interaction [6, 7] and hence the present calculation is consistent with those results. 

III. INTERMEDIATE irN AMPLITUDE 

The combination of Figs. 2 and 3 gives rise to the TPE-3NP, associated with the six 
diagrams shown in Fig. 4. In the sequence, we discuss their individual contributions to the 
subamplitudes D and B . We are interested only in the longest possible component of 
the potential and numerators of expressions are systematically simplified by using k 2 — > fi 2 
and k' 2 — > /i 2 . In configuration space, this corresponds to keeping only those terms which 
contain two Yukawa functions and neglecting interactions associated with Figs. 1 (c) and 1 
(d). 

• diagrams (a) and (6) : The crosses in the nucleon propagators of Figs. 4 (a) and 4 (b) 
indicate that they do not include forward propagating components, so as to avoid double 
counting when the potential is used in the dynamical equation. The covariant evaluation of 
these contributions is based on Eqs. (10-13). Denoting by p the momenta of the propagating 
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(a) 



(b) 



(c) 



» « + 
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(d) 



(e) 



(f) 



FIG. 4: (Color online) Structure of the 0(q 4 ) two-pion exchange three-nucleon potential 



nucleons, the factors l/(s — m 2 ) and l/(u — m 2 ) are decomposed as 

1 1 1 



(20) 



(p°) 2 -E 2 2E(p°-E) 2E(p° + E) ' 

with E = \Jm 2 +p 2 . The first term represents forward propagating nucleons, associated with 
the iteration of the OPEP, whereas the second one gives rise to connected contributions. 
Discarding the former and using the results of Appendix A, one has 

V({«} - m 2 ) -> -V [4™ 2 + {3q 2 + q 2 p /3 + 16Q 2 j3 ± l0q r -Q p /V3T2q p -Q r /Vs)] ■ (21) 



After appropriate truncation, one obtains 



D 



B 



ab 



D 



0{q 



rrr 

2\ 



(2f-t) - 0{q 2 ) 



ab 



D 



ab 



2m 2 
0(q 2 ) , 



v - 0{q 2 ) 



(22) 
(23) 
(24) 
(25) 



where we have used the fact that, in the case of virtual pions, v —± 0(q 2 ). 

• diagrams (c) and (d): These contributions are purely polynomial, can be read directly 



9 



from Eqs. (16-19), and are given by 



4 Cl 



P 



c 3 



9aM 



(2/x 2 -t)^0(g 2 



b: 



1 



cd 



v -> 0(^ 



1 + 2c 4 m _ ^m/x _^ 



(26) 
(27) 
(28) 

(29) 



cd 2/2 ' / 2 8tt/4 
diagrams (e) and (/): The medium range components of the intermediate ttN amplitude 



are 



d: 



D 



b: 



9aV 
64tt 2 /* 

> 0(q A ) 



g A m\i 



where lit is the dimensionless Feynman integral 
fi 2 F(a) 



(2t-/i 2 ) [(l-t/2/i 2 )n 4 -27r] ^0(q 3 ) , 



[(l-t/V)II t -7r] -0(g) 



n, 



da 



t-M 2 



M = 2/i/a , F(a) = — tan" 1 



ma 



L/*(l-oV2) 



(30) 
(31) 
(32) 

(33) 



The amplitude proportional to z/, is 0(q 3 ) for free pions and here becomes 0(q A ). Thus, 
in fact, diagram (/) does not contribute to the TPE-3NP at 0(q A ). 

• full results: The Golberger-Treiman relation g/m = qaI fn is valid up to 0(q 2 ) and can 
be used in diagrams (a) and (b). One then has 

9a( 1 +9a)^ 9a» 



D 

where 



a(2/i 2 ) , (2/i 2 -t) f g\ 



J 7T 



+ 



P 



5m 



C 3 + 



a(t = 2/i 2 



16tt/ 2 
-4 ci /i 



1287T 2 / 2 



;i-2*/^)n t 



2 _ 3fl> 3 
32tt/ 2 



(34) 



(35) 



is the value of the scalar form factor at the Cheng-Dashen point [14]. The remaining ampli- 
tudes read 



5 J 
D~ 

B~ 



v 



0(q 2 ) 
1-9 A 

VI 

1 + 4 c^m 
VI 



167T/4 



9 A m ^ 

16 tt 2 /* 



;i-t/4/i 2 )n t 



(36) 
(37) 

(38) 



10 



The subamplitudes D ± and B ± begin at 0(q 2 ) and one needs just the leading terms in the 
spinor matrix elements of Eq. (8), which is rewritten as 



T+ = 2m D' 



T~ = 2m D- + KT {3) -k'xkB- 



(39) 
(40) 



with D + -> 0(q 2 ) + 0(q 3 ), D~ -> 0(g 2 ), and £" -> C(g°)+C>(g). 

• C?(<7 3 ) reduction: In order to compare our amplitudes with previous 0(q 3 ) results, one 
notes that, in case corrections are dropped, one would have 

<r(0) (2/i 2 -t) 



D" 



Br 



f 2 


J TV 


1 


2 c^m 


Uf 2 _ 


f2 
J n 



8m 



c 3 



(41) 
(42) 



These expressions agree with those derived directly from a chiral lagrangian [24], except 
for the terms within square brackets in both D + and B~ . The former corresponds to a 
Born contribution whereas the latter is due to diagram (c) in Fig. 4, associated with the 
Weinberg- Tomozawa term. 



IV. TWO-PION EXCHANGE POTENTIAL 

The expansion of the TPE-3NP up to 0(q A ) requires only leading terms in vertices and 
propagators. In order to derive the non-relativistic potential in momentum space, one divides 
Eq. (7) by the relativistic normalization factor \J2E ~ \/2m for each external nucleon leg 
and writes 1 



4/2 k 2 +f i 2 k' 2 +^ 



T W. T (^+_ iT (D XT (2). T (3)( r + _L (r (3). ji ./ xfcF 

2m 



The configuration space potential has the form 



(43) 



V 3 (r, p) = r (1) -r (2) V 3 + (r, p) + r (1) x r (2) -r (3) Vf(r, p) + cyclic permutations, (44) 



1 One notes that this expression is identical with Eq. (33) of Ref. [10] divided by 8m 3 
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with 



V+(r,p) = C+ ct^-x 31( t^-x 23 U x {x 3l )U x {x 23 ) 

+ C+ {(l/9)o-W-o-( 2 ) [U(x 31 )-U 2 (x 31 )\ [U(x 23 )-U 2 (x 23 )\ 
+ (1/3)<tW-x 23 <tW-x 23 [U(x 31 )-U 2 (x 31 )] U 2 (x 23 ) 
+ (1/3) ctW-x 31 <tM -x 31 U 2 (x 31 ) [U{x 23 )-U 2 {x 23 )\ 
+ (T {1) -x 31 cr (2) -a;23 x 31 -x 23 U 2 (x 31 ) U 2 (x 23 )} 

+ C+ a W • <r« • V 2 ' 3 Vi ■ V^ 3 [/° - 2 J 1 ] , (45) 

V 3 -(r, P ) =Cr{ (1/9) ^W^-O^ [f/(x 31 )-[/ 2 (x 31 )] [*7(Z23)-EMZ23)] 

+ (l/3)^ 3 W 1 )-£ 23CT ( 2 )-£ 23 [C/(x 3 i) -£/ 2 (x 31 )] f/(x 23 ) 
+ (l/3)<r«-& 31 <r( 2 W 3 )-& 31 E/ 2 (x 3 i) [C/(x 23 )-^ 2 (x 23 )] 
+ <T (1) -i 3 i <r (2) -a; 2 3 (T (3) i 3 iXX23 f/ 2 (x 3 i)t/ 2 (x 23 )} 
+ C 2 -{ -W-(iV 3 f-iV 2 f) ^-i 23 [U(x 31 )-U 2 (x 31 )} Ui(x*) 
+ tr^-x ja tr^-(iVg-iVg) C^si) [t/(x 23 )-[/ 2 (x 23 )] 

+ 3<r (1) - : z 31 <7 (2) -£ 23 (W 3 f -iV 2 f ) • [x 31 U 2 (x 31 ) f/i(x 23 ) + Ui(x 31 ) U 2 (x 23 )]} 
+ C 3 " tr W • V/, <r® . v i CT (3) . v / x v i [7 o _ 7 i /4] . (46) 



The profile functions are written in terms of the dimensionless variables x^ = fi and read 



U(x) 



x 



f 1\ 



f/ 2 (x) 



x/ a; 
3 3\ e" 1 



X X* I X 



^ 2 J (2tt) 3 (2vr) 3 



fc 2 + u 2 fc /2 + u 2 



(47) 
(48) 

(49) 

n t (t) . (50) 
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The last function involves the loop integral given in Eq. (33) and is discussed further in 
Appendix C. The gradients act on the functions I n , whereas the V™- act only on the 
wave function and give rise to non-local interactions, as discussed in Appendix D. 

The strength coefficients are the following combinations of the basic coupling constants 

C ' + = 6W ff(2 " 2) ' (51) 

a = (jA + mC3 + Mgpt ) , (52) 

32vr 2 /4 m y 8 16tt/ 2 / 

C * = 4096 7T3/6 ' (53) 

<t TV 

256 7r 2 f^m \ 8tt/ 2 J 

Ci = • (55) 



768 7r 2 /4 m 
2048 tt 3 /6 



^ = "^#4^- (56) 



V. STRENGTH COEFFICIENTS 

The strength constants of the potential involve a blend of four well determined param- 
eters, namely m = 938.28 MeV, n = 139.57 MeV, g A = 1.267 and f n = 92 A MeV, with 
the scalar form factor at the Cheng-Dashen point and the LECs C3 and C4, which are less 
precise. As far as <r(2/i 2 ) is concerned, we rely on the results [25] <r(2/z 2 ) — <t(0) = 15.2 ± 0.4 
MeV, o-(O) = 45 ±8 MeV, and adopt the central value <r(2/i 2 ) = 60 MeV. The values quoted 
for the LECs in the literature vary considerably, depending on the empirical input employed 
and the chiral order one is working at. A sample of values is given in Table I. 

Our work is based on the 0(q 3 ) expansion of the intermediate irN amplitude and, for 
the sake of consistency, one must use LECs extracted at the same order. The kinematical 
conditions of the three-body interaction are such that the variable v is 0(q 2 ), an order of 
magnitude smaller than the threshold value, v — fi. This makes information encompassed in 
the subthreshold coefficients better suited to this problem and we use results from Appendix 
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TABLE I: Some values of the LECs C3 and C4; m is the nucleoli mass. 



Reference Chiral order irN input m C3 m C4 





3 


amplitude at v — 0, t = 


-5.00 ± 1.43 3.62 ±0.04 


[26] 


3 


amplitude at v = 0, t = 2/i 2 /3 


-5.01 ± 1.01 3.62 ±0.04 


[27] 


3 


scattering amplitude 


-5.69 ±0.04 3.03 ±0.16 


[15] 


4 


subthreshold coefficients 


-3.4 


2.0 


[15] 


4 


scattering lengths 


-4.2 


2.3 


tree 


2 


subthreshold coefficients 


-3.6 


2.0 


this work 


3 


subthreshold coefficients 


-4.9 


3.3 



B in order to write 

,2r+ 9A m V 77 9A m V /k«N 

mc 3 = -mUd m - — - -^-^ , (57) 

mC4_ 2 4 + 16 tt/? ' { } 

Adopting the values for the subthreshold coefficients given by Holder [13], namely dji = 
1.14 ± 0.02 /i -3 and 6q = 10.36 ± 0.10 fi~ 2 , one finds the figures shown in the last row of 
Table I. These, in turn, produce the strength coefficients displayed in Table II. For the sake 
of comparison, we also quote values employed in our earlier calculation [10] and in two TM' 
versions [28] of the Tucson- Melbourne potential [8]. 

TABLE II: Strength coefficients in MeV. 



reference 


ct 


ct 


^1— 

°3 °1 °2 


C3 


this work 


0.794 


-2.118 0.034 0.691 0.014 


-0.067 


Brazil [10] 


0.92 


-1.99 


- 0.67 - 




TM'(93) [28] 


0.60 


-2.05 


- 0.58 - 




TM'(99) [28] 


0.91 


-2.26 


- 0.61 - 





Changes in these parameters represent theoretical progress achieved over more than two 
decades and it is worth investigating their origins in some detail. With this purpose in 
mind, we compare present results with those of our previous 0(q 3 ) calculation [10]. At the 
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chiral order one is working here, new qualitative effects begin to show up, associated with 
both loops and non-local interactions. They are represented by terms proportional to the 
coefficients C% , and C% in Eqs. (45) and (46). 

The irN coupling is now described by g\^ 2 / f 2 = 3.66 whereas, previously, the factor 
g 2 fi 2 /m 2 = 3.97 was used. From a conceptual point of view, the latter should be preferred, 
since g is indeed the proper coupling observable. In chiral perturbation theory, the difference 
between both forms is ascribed to the parameter Act = — 2disfi 2 /g, which describes the 
Goldberger-Treiman discrepancy [15]. As this is a 0(q 2 ) effect, both forms of the coupling 
become equivalent in the present calculation. On the other hand, the empirical value of g is 
subject to larger uncertainties and the form based on is more precise. Our present choice 
accounts for a decrease of 8% in all parameters. 

The relations C s , C p and C*f «-> —C allow one to compare Eqs. (45) 

and (46) with Eq. (67) of Ref. [10]. One notes that the latter contains an unfortunate 
misprint in the sign of the term proportional to C , as pointed out in Ref. [29]. In the 
earlier calculation, the coefficient C s was based on a parameter [30] a a = 1.05/i -1 , which 
corresponds to <r(2/i 2 ) = 64 MeV. The results of Table II show that the values of and 
C-f are rather close to those of C v and —C' p . This can be understood by rewriting Eqs. (52) 
and (54) in terms of the subthreshold coefficient dji an d &oo as follows 

(59) 
(60) 

Numerically, this amounts to C 2 + = -(1.845 + 0.110 + [0.163]) MeV and C7/ = (0.624 + 
[0.067]) MeV. The second term in the former equation was overlooked in Ref. [10] and 
should have been considered there. The square brackets 2 correspond to next-to-leading 
order contributions and yield corrections of about 8% and 11% to the leading terms in 
and C-f, respectively. 3 As the model used in Ref. [10] was explicitly designed to reproduce 
the subthreshold coefficients quoted by Hohler [13], it produces the very same contributions 
as the first terms in Eqs. (59) and (60). 

2 These factors can be traced back to loop diagrams in Fig. 3 and are dynamically related with the term 
proportional to C^, as we discuss in Appendix C. 

3 When comparing the new coefficients with those in the second row of Table II, one should also take into 
account the 8% effect due to the Goldberger-Treiman discrepancy. 
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VI. NUMERICAL RESULTS FOR THREE-NUCLEON SYSTEMS 



In order to test the effects of the TPE-3NP at 0(q A ), in this section, we present some 
numerical results of Faddeev calculations for three-nucleon bound and scattering states. The 
calculations are based on a configuration space approach, in which we solve the Faddeev 
integral equations [31-33], 

1 



"12,3 



E + ie — H — V 12 
x [V 12 ($i + $ 2 ) + W 3 ($i + $ 2 + $ 3 )] , 

(and cyclic permutations), (61) 

where H 12j 3, which does not appear in the bound state problem, is an initial state wave 
function for the scattering problem, H is a three-body kinetic operator in the center of 
mass, V\2 is a nucleon-nucleon (2NP) potential between nucleons 1 and 2, and W% is the 
3NP displayed in Fig. 2. Partial wave states of a 3N system, in which both NN and 3N 
forces act, are restricted to those with total NN angular momenta j < 6 for bound state 
calculations, and j < 3 for scattering state calculations. The total 3N angular momentum 
(J) is truncated at J = 19/2, while 3NP is switched off for 3N states with J > 9/2 for 
scattering calculations. These truncation procedures are confirmed to give converged results 
for the purposes of the present work. 

When just local terms are retained, i% in Eq. (43) can be cast in the conventional form 
[8-10] 

x "(r (1) •r (2) ){a + 6(fc-fc / )} 

-(zr (1) x r (2) • r (3) ) (ia® ■ k' x k)d I , (62) 

where the coefficients, a, b, and d are related with our potential strength coefficients by 

[C+ C 2 + , Cr] = j^y 2 [-a/i 4 , 6/i 6 , -rf/i 6 ] . (63) 

The values of the coefficients, a, b, and d for the TPE-3NP at 0(q A ) are shown in Table III, 
as BR-(9(g 4 ). In this table, the values for the older version of the Brazil TPE-3NP, BR(83) 
[10], and the potential up to 0(q 3 ) given by Eqs. (41-42), BR-C(g 3 ), are shown as well. 
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TABLE III: Coefficients a, b, and d of the TPE-3NP. 



3NP 


a fj, 


b fi 3 


d fjP 


BR-0(g 4 ) 


-0.981 


-2.617 


-0.854 


BR-0(g 3 ) 


-0.736 


-3.483 


-1.204 


BR(83) 


-1.05 


-2.29 


-0.768 



In Eq. (62), the function F(k 2 ) represents a 7rNN form factor. We apply a dipole form 
factor with the cut off mass A, > which modifies the profile functions U(x), U\(x), 

and C/ 2 (a;) in Eqs. (47-49) as 




with A = A/fi. 

We choose the Argonne Vi§ model (AV18) [34] for a realistic NN potential, by which the 
triton binding energy (B3) becomes 7.626 MeV, underbinding it by about 0.9 MeV compared 
to the empirical value, 8.482 MeV. As it is well known, the introduction of the TPE-3NP 
remedies this deficiency. The amount of attractive contribution depends on the cutoff mass 
A, as shown in Fig. 5. The solid curve shows the dependence of B 3 on A for the calculation 
with the BR-C(g 4 ) 3NP in addition to the AV18 2NP (AV18+BR-C(g 4 )). In the figure, 
the empirical value and the AV18 result are displayed by the dashed and dotted horizontal 
lines, respectively. Due to the strong attractive character of the 3NP, B3 is reproduced by 
choosing a rather small value of A, namely 660 MeV. In the same figure, the A-dependence 
of B 3 for AV18+BR-0(g 3 ) is displayed by a dashed curve and that for the AV18+BR(83) by 
a dotted curve. From these curves we see that AV18+BR-(9(g 3 ) reproduces B 3 for A = 620 
MeV and AV18+BR(83) for A = 680 MeV. In other words, the BR-0(g 4 ) 3NP is slightly 
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more attractive than the BR(83) 3NP and a large attractive effect occurs when one moves 
from the TPE 0(q^) 3NP to the 0(q 3 ) 3NP. This tendency is strongly correlated with the 
magnitude of the coefficient b, as shown in Table III. This can be understood as a dominant 
contribution to B 3 from the component of the TPE-3NP associated with the coefficients b. 
This dominance is shown in Table IV, where we tabulate calculated B 3 for the AVI 8 plus 
the BR-0(g 4 ) 3NP and plus each term of the BR-0(g 4 ) coming from the coefficients a, b, 
and d. 




FIG. 5: (Color online) The triton binding energy B3 as functions of the cutoff mass A of the irNN 
dipole form factor. The solid curve denotes the result for AV18+BR-0(g 4 ), the dashed curve for 
AV18+BR-C>(g 3 ), and the dotted curve for AV18+BR(83). The horizontal lines denote the AV18 
result (dotted line) and the empirical value (dashed line). 

In Fig. 6, we compare six calculated observables for proton-deuteron elastic scattering, 
namely differential cross sections cr(9), vector analyzing powers of the proton A y (9) and 
of the deuteron iTu(9), and tensor analyzing powers of the deuteron T 20 (6 I ), T 2 i(9), and 
^22(0), at incident proton energy E l § b = 3.0 MeV, (or incident deuteron energy E l f b = 6.0 
MeV,) with experimental data of Ref. [35, 36]. In the figure, the solid curves designate 
the AV18 calculations and the dashed curves the AV18+BR-(9(g 4 ) calculations, which are 
almost indistinguishable from the AV18+BR-(9(g 3 ) and AV18+BR(83) calculations, once 
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TABLE IV: Triton binding energy for the AV18 2NP plus the BR-C(g 4 ) 3NP for each term of the 
BR-C(g 4 ) 3NP with A = 660 MeV. AB3 means the difference of the calculated binding energy 
from that of the AV18 calculation. 



B 3 (MeV) AB 3 (MeV) 



AV18+BR-C(g 4 ) 8.492 0.866 

AV18+BR-0(g 4 )-a 7.673 0.047 

AV18+BR-C(g 4 )-6 8.241 0.615 

AV18+BR-0(g 4 )-d 7.787 0.161 



the cut off masses are chosen so that B 3 is reproduced. 

It is reminded that the TPE-3NF gives minor effects on the vector analyzing powers. 
This happens because the exchange of pions gives essentially scalar and tensor components 
of nuclear interaction in spin space, which are not so effective to the vector analyzing powers. 
On the other hand, as is noticed in Refs. [37, 38], at E l £ b = 3.0 MeV, the TPE-3NP gives a 
wrong contribution to the tensor analyzing power T 2 \{9) around 9 = 90°. 

In Fig. 7, we compare calculations of observables in neutron-deuteron elastic scattering 
at E l $ b = 28.0 MeV with experimental data of proton-deuteron scattering Ref. [39]. At 
this energy, discrepancies between the calculations and the experimental data in the vector 
analyzing power iT n (9) appear at 9 ~ 100°, where iT n (9) has a minimum, and at 9 ~ 140°, 
where %T\\ (9) has a maximum, which are not compensated by the introduction of the TPE- 
3NP. On the other hand, while the AV18 calculation almost reproduces the experimental 
data of T2\(9) at 9 ~ 90°, the introduction of the TPE-3NP gives a wrong effect, as in the 
Eft 3 = 3 MeV case. 

These results set the stage for the introduction of terms associated with the coefficients 
C%, C2 1 and C3", Eqs. (44-45), which are new features of the 0{q i ) expansion of the TPE- 
3NP. Terms proportional to C^, which include the rather complicated function 1(7*31,^3) 
given in Appendix C, arise from a loop integral, Eq. (33). On the other hand, the term with 
C2 corresponds to a non-local potential and includes the gradient operator V*" , which acts 
on the wave function and arises from the kinematical variable v. Both kinds of contributions 
are not expressed in the conventional local form shown in Eq. (62), which involves only 
the coefficients C^~, C£, and C-f, and the full evaluation of their effects would require an 
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E = 3.0 MeV (E = 6.0 MeV) 

N 1 d * 




60 120 180 60 120 180 60 120 

8 (deg) 6 (deg) 8 (deg) 



FIG. 6: (Color online) Proton-deuteron elastic scattering observables at E l ^ b = 3.0 MeV. Solid 
curves are calculations for the AV18 potential, and dashed curves for the AV18+BR-C(g 4 ). Ex- 
perimental data are taken from Refs. [35, 36]. 

extensive rebuilding of large numerical codes. However, the coefficients of the new terms are 
small, and in this exploratory paper we estimate their influence over observables as follows. 

The function I(r 31 , r 23 ) is approximated by Eq. (Cll), which amounts to replacing IL(t) 
by a factor — ir. Further, the kinematical factors in front of Ht(t) in Eqs. (34) and (38), 
namely 1 — 2t/fi 2 and 1 — t/4fi 2 , are approximately evaluated by putting t m 2/i 2 , which 
yields —3 and 1/2, respectively. By this procedure, the coefficients C% and C% are absorbed 
into C£ and Cf, or in b and d respectively, and one has 

AC+ = -3C+, ACT = 072- (67) 

Numerically, this corresponds to AC 2 + = -0.102 MeV ~ ^C 2 + and ACf = -0.034 MeV ~ 
— ^jCr , or Ab = — 0.125(/x -3 ) and Ad = 0.042(/i -3 ). The net change produced in the triton 
binding energy is +0.026 MeV (+0.037 MeV from AC 2 + and -0.011 MeV from ACf), just 
about 1/30 of the total increase in B3 due to the local terms of the BR-(9(g 4 ) TPE-3NP. 

The non-local term proportional to is more involved and we restrict ourselves to a 
rough assessment of its role. We replace the variable v by a constant (y) and assume, for 
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E J ab = 28.0 MeV (E ' ab = 56.0 MeV) 

N d 




9 cm (deg) 9 cm (deg) 9 cm (deg) 

FIG. 7: (Color online) Nucleon-deuteron elastic scattering observables at E l $ b = 28.0 MeV. Curves 
are calculations for neutron-deuteron scattering. Solid curves denote calculations for the AV18 
potential and dashed curves for the AV18+BR-C(g 4 ). Experimental data are those for proton- 
deuteron scattering taken from Ref. [39]. 

2 

example, that (u) = j-. This changes the term in Eq. (46) into the very simple form 
Vf{r, p) = Cr (• • • ) + i(7 2 -<7« • x 31 <t® . x 23 U 1 (x 31 )U 1 (x 23 ) + C 3 (• • • ) , (68) 

with 

C- = - Al^A {v ) = = 0.021 MeV . (69) 

2 4/2 2/2 V '{Anf 512tt2/> ^ } 

Except for the isospin factor, this term is similar to that with (or a), which adds 
about 0.05 MeV to the triton binding energy. Since the potential strength C 2 is about 3 % 
of Ci , its contribution to the binding energy may be estimated to be a tiny 0.001 MeV. 

VII. CONCLUSIONS 

In the framework of chiral perturbation theory, three-nucleon forces begin at 0(q 3 ), with 
a long range component which is due to the exchanges of two pions and relatively simple. 
At 0(q 4 ), on the other hand, a large number of different processes intervene and a full 
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description becomes rather complex. For this reason, here we concentrate on a subset of 
0(q A ) interactions, namely that which still involves the exchanges of just two pions. This 
part of the 3NP is closely related with the 7rN amplitude, and the expansion of the former 
up to 0(q A ) depends on the latter at 0(q 3 ). 

Our expressions for the potential are given in Eqs. (44-56) and the new chiral layer of the 
TPE-3NP considered in this work gives rise to both numerical corrections to strength coef- 
ficients of already existing terms (C± , C£, Cj~) and new structures in the profile functions. 
Changes in numerical coefficients lay in the neighborhood of 10% and can be read in Tables 
II and III. New structures, on the other hand, arise either from loop functions representing 
form factors or the non-local terms associated with gradients acting on the wave function. 
They correspond to the terms proportional to the parameters and C^, which are 

small and compatible with perturbative effects. 

In order to insert our results into a broader picture, in Table V we show the orders at 
which the various effects begin to appear, including the drift potential derived recently [40] . 



TABLE V: Chiral picture for two- and three-body forces. 



beginning 


TWO-BODY 


TWO-BODY 




THREE-BODY 


O(q ) 


0PEP:y T -,C 5 - s 








0{q 2 ) 


OPEP: V D 


TPEP: V^;V+, 






0{q 3 ) 




TPEP:E- s ,Ef 


v~ ■ v + v + 


TPEP: Ci;C+,C£ 


0(q 4 ) 




TPEP: VjjJ; Vq . 




TPEP: C^; C%, C% 



The influence of the new TPE-3NP over three-body observables has been assessed in both 
static and scattering environments, adopting the Argonne Vis potential for the two-body 
interaction. In order to reproduce the empirical triton binding energy, the 0(q 4 ) potential 
requires a cutoff mass of 660 MeV. Comparing this with the value of 680 MeV for the 1983 
Brazil TPE-3NP, one learns that the later version is more attractive. 

In the study of proton-deuteron elastic scattering, we have calculated cross sections cr(9), 
vector analyzing powers A y {9) of the proton and 2X11(6?) of the deuteron, and tensor analyzing 
powers T 2 o(9), T 2 i(6>), and T 2 2{9) of the deuteron, at energies of 3 and 28 MeV. Results are 
displayed in Figs. 6 and 7, where it is possible to see that there is little sensitivity to 
the changes induced in the strength parameters when one goes from 0(q 3 ) to 0(q 4 ). Old 
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problems, as the A y (9) puzzle, remain unsolved. 

The present version of the TPE-3NP contains new structures, associated with loop inte- 
grals an non-local operators. Their influence over observables has been estimated and found 
to be at least one order of magnitude smaller than other three-body effects. A more detailed 
study of this part of the force is being carried on. 

APPENDIX A: KINEMATICS 

The coordinate describing the position of nucleon i is and one uses the combinations 
R= (r 1 + r 2 + r 3 )/3, r = r 2 -r 1: p = (2 r 3 -r 1 -r 2 )/ v / 3 , (Al) 

which yield 

The momentum of nucleon % is p { and one defines 

P = Pi+P 2 +P 3 , Pr = (P 2 -Pi)/2, P p = (2p 3 -p 1 -p 2 )/2 v / 3. (A3) 
Initial momenta p and final momenta p' are used in the combinations 

Q = (P'+P)/2, q = (P'-P), (A4) 

Qr = (Pr+Pr)/ 2 , Qr = (Pr~Pr) > ( A5 ) 

Q P = (p' p +P p )/2 , % = (P P ~P P ) ■ (A6) 
In the CM, one has P = and the three-momenta are given by 

Pi = -{Q r -q r /2) - (Q p -q p /2)/V3 , p[ = -(Q r +q r /2) - (Q p + q p /2)/v / 3 , (A7) 
p 2 = (Q r -q r /2) - (Q p -q p /2)/V3 , p' 2 = (Q r + q r /2) - (Q p +q p /2)/V3 , (A8) 
p 3 = 2(Q p -q p /2)/V3 , pf 3 = 2(Q p + q p /2)/V3 . (A9) 

Energy conservation for on-shell particles yield the non-relativistic constraint 

Q r -q r + Q P -Q P = o. (Aio) 
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The momenta of the exchanged pions are written as 

k=p 1 -p' 1 , k'=p' 2 -p 2 , (All) 
k° = -(q r + q p /V3)-(Q r + Q p /V3)/m , k = q r + q p /V3 , (A12) 
k'° = (q r -q p /V3)-(Q r -Q p /V3)/m , k' = q r -q p /V3 , (A13) 

and the Mandelstam variables for nucleon 3 read 

8 = (p 3 + k) 2 = m 2 - (q r + q p /V3) ■ (q r + 2 Q r -g p /v / 3 + 2 V / 3 Q p ) + 0{q A ) , (A14) 
u = (p 3 -k') 2 = m 2 - (q r -q p /V3) ■ (q r + 2 Q r + q p /v / 3-2 V / 3 Q p ) + 0(q 4 ) , (A15) 
v = ( s -u)/4m = -2q r -Q p /V3 + 0(q 4 ) . (A16) 

In the evaluation of the intermediate amplitude, one needs 

[u(p') w(f>)] (3) ~ 2m + 0(q 2 ) , (A17) 
[^-u(jpf) a^(j/-pYK v u(p)]^ ~2i<r®-q p xq r /V3 + 0(q 4 ) . (A18) 

The ttN vertex for nucleon 1 is associated with 

[n(p')7 5 n(p)] (1) ~ crU-(q r +q p /V3) + 0(q 3 ) , (A19) 

and results for nucleon 2 are obtained by making q r — > —q r . 

APPENDIX B: SUBTHRESHOLD COEFFICIENTS 

The polynomial parts of the amplitudes T R , Eqs. (30-35), are determined by the sub- 
threshold coefficients of Ref. [15]. The terms relevant to the 0(q 3 ) expansion are written as 
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where the parameters q and di are the usual coupling constants of the chiral lagrangians 
of order 2 and 3 respectively [41] and the tilde over the latter indicates that they were 
renormalized [15]. Terms within square brackets labeled (mr) in these results are due to 
the medium range diagrams shown in Fig. 3 and have been included explicitly into the 
functions D^ r and B^ r . Terms bearing the (WT) label were also explicitly considered in 
Eqs. (15-19). The subthreshold coefficients are determined from ttN scatterig data and a 
set of experimental values is given in Ref. [13]. 



APPENDIX C: FUNCTIONS I n 



The functions I n , describing loop contributions, are given by 



I n (r 3 i,r 23 ) 



167T 



dk dk Mk-r 3 i+k'-r 2 3) 
(2tt) 3 (2tt) 3 



t 

7 2 



k 2 +fi 2 k' 2 +fi 2 



Using the definition Eq. (33) and the Jacobi variables Eq. (Al), one writes 



/ n (r 31 ,r 23 ) 



4v; 

Tp? 



-f( r 3i' r 23) = 1287T / da tan 1 
Jo 



ma — a? 



L(a;r,p) 



/i(l-a 2 /2) 
dq dQ e <Q-r-V3q-p/2) 



L(a;r,p) 
1 



1 



n t (t) . (ci) 

(C2) 
(C3) 
(C4) 



(2tt) 3 (2vr) 3 aV + 4/i 2 [(Q-qf+y?] [{Q +q )2 + ^ • 
The numerical evaluation of the function L is can be simplified by using alternative repre- 
sentations. 
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form 1: One uses the Feynman procedure for manipulating denominators, which yields 

f dq dQ e i(Q-r-Vsq-p/2) 1 



L(a;r,p) = db J 



(2tt) 3 (2tt) 3 a 2 q 2 + V [(Q 2 + q 2 /4+/i 2 )-(l-2%-Q] 



2 



If 1 f dq e *[(i-2&)r-V3p]-g/2 e -e, 
d6 



8 Wo J (2tt) 3 a 2 g 2 +V 



6 = ^2+6(1-6) q 2 . (C5) 
Performing the angular integration over q, one has 

Tt \ 1 f\ h I A e ' 0r smq[(l-2b)r-V3p]/2 

Lla; r, p) = / do dq q ^ . _ . — — — = . (Co) 

Kl ^' 16^J J qq e (a 2 qr 2 +V) [(l-26)r- >/3p]/2 1 ; 

• form 2: The Fourier transform 

= / daj e-*- 25 ^ (C7) 

allows one to write 

1 1 C P -V\ r 3i+Z\ -fi\r 23 -Z\ p -2nz/a 

L(a;r,p) = — - dz \ . (C8) 

These results may be further simplified by means of approximations. 

• heavy baryon approximation: In the limit m — > oo, corresponding to the heavy 
baryon case, one uses F(a) — > Ait/a 2 in Eq. (33) and Eqs. (C5) and (C7) yield, respectively, 



^(^3i, ^23) ~ - db dq 



2 f\, f°°, , q 1 e~ @r sing [(1-26) r-V?>p}/2 

tan 



7T 



JO 



2/' 

I r e -Ai|r 3 i+z| e -v\r 23 -z\ e -2 M 

- dz j j ; 

J |T*31+^| \r23-z\ 2flZ 



p6 [(l-26)r-V3p]/2 



(C9) 



J(r 31 , r 23 ) ~ - j dz ,_ ; _, — ^-^ . (CIO) 



• multipole approximation: The integrand in Eq. (CIO) is peaked around z = and a 
multipole expansion of the Yukawa functions produces 

/(r-si, r 23 ) ~ U(x 31 ) U(x 23 ) + ■■■ . (Cll) 

The same result can also be obtained by using the expansion Tlt{t) ~ — 7r[l + t/12/i 2 + 
t 2 /80/i 4 + •••], valid for low t, directly into Eq. (CI). 
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APPENDIX D: NON-LOCAL TERM 



In configuration space, the variable Q p corresponds to a non-local operator, represented 
by a gradient acting on the wave function. In order to make the dependence of £3 on Q p 
explicit, one writes 

h = [Q P ]iXi(q r ,q p ) , (Dl) 
where X is a generic three-vector, and evaluates the matrix element 



1 




[(27T). 





X 



e i[QAr>~r) + Q^-p )+ q A r> + r)/2 +qp . { p>+p V 2] ^ Q ^ q %) 
1 



dr dp 



-V p r(r,p) 



.(2t)J 

x | dg p dg p e 4 [9r-r+g p .p] X,(q r ,q p ). 
This yields the potential 

[2/V3] 3 



-V p tp(r,p) 



Vs(r,p) 



(2*) 



(D2) 



(D3) 



where the operator V = V — V acts on/?/ on the wave function. An alternative form can be 
obtained by integrating Eq. (D2) by parts, and one finds 

[2/V3] 3 



V 9 (r,p) = - 



(2tt)< 



J dq r dq p A«r-r+QM X (q r , q p ) • [- i V/ ] 



V p • dq r dq p A<lr- r +<lM X(q r ,q^ 



(D4) 



In the case of the three-body force, the only non-local contribution is associated with the 
subamplitude D~, Eq. (37), which yields 



k z +fi 2 k' 2 +fi 



9 2 a(9 2 a-1) 



(k'+k)i . (D5) 



The action of V p on the second term of Eq. (D4) gives rise to an integrand proportional to 
(k 2 — k 2 ), which has short range and does not contribute to the TPE-3NP. Therefore it is 
neglected. 
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